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Abstract. We examine Pr 3+ crystal field models with near-degeneracy of the two 
lowest crystal field levels, where interaction-induced quadrupolar and dipolar moments 
are as important as the permanent moments of the crystal field ground state. We find 
that the Tx~ T4 level scheme yields a successful description of the antiferroquadrupolar 
ordering of PrFe4Pi2- For PrBa2Cu306, we argue that quadrupolar interaction is 
important for understanding the Pr ordering transition at 11K. 



1. Introduction 

The highly degenerate /-shells of rare earth and actinide ions support a great variety 
of local degrees of freedom: magnetic (dipolar), quadrupolar, and octupolar. For 
localized /-electron systems, these (and some of their combinations) are potential order 
parameters, leading to rich phase diagrams, as for CeB 6 p. Field-induced octupoles are 
known to influence the shape of the phase boundaries for CeB 6 or PrPb3 [12] , but it is 
an outstanding question whether breaking time-reversal invariance with purely octupole 
order is possible; recent observations on NpC>2 indicate that it may be so 0. The /-shell 
multipolar degrees of freedom may also interact and mix with those of other electrons: 
the unusual heavy fermion behavior of PrFe4Pi2 0, an d the exotic superconducting 
phases of PrOs 4 Sbi 2 [H] have been tentatively ascribed to Pr 3+ quadrupole fluctuations. 

Detailed understanding has been achieved for CeBg and some other Ce-based 
systems. The reason is that Ce 3+ has only a single /-electron, and therefore microscopic 
analysis is here the easiest. Of particular interest is cubic CeBg in which the fourfold 
degenerate r 8 ground state supports the full array of 3 dipoles, 5 quadrupoles, and 7 
octupoles pp. A microscopic analysis P indicates that octupole-octupole interactions 
should not be weaker than quadrupole-quadrupole interactions, which gives rise to a 
close competition between different kinds of order. In fact, if all interactions were exactly 
equal, the model would be SU(4)-symmetrical, which sparked early theoretical interest 
in the subject [4... The realistic multipolar model is not fully symmetrical, but the 
nearness of the SU(4) point is still felt in fluctuations [7^ Whether the ground state of 
an SU(4) model is long range ordered, or it is rather liquid-like (the SU(4) version of 
RVB), depends on details like the lattice structure, and the range of the interaction [Hj. 



Quadrupolar interactions in Pr compounds: PrFe^Pyi and PrBa^ Cu% 0§ 



2 



The approximate realization of high symmetries like SU(4) is of current interest for d- 
and /-electron systems alike. 

We expect that the result that various multipolar couplings are of the same order 
of magnitude, holds generally for light rare earths. The question arises what can we say 
about 4/ 2 systems. 

Pr compounds show varied behavior. A few examples: PrPb3 undergoes an 
antiferroquadrupolar transition but it does not order magnetically, though both dipolar 
and octupolar interactions are thought to be important ^2]- PrBg has two magnetic 
phase transitions but the magnetic anisotropy is explained by quadrupolar interactions 
[TT] . PrBa2Cu307_a; is the only member of the group (RE)Ba2Cu307 (RE=rare earth), 
which is not a high-T^ superconductor [TH]. It is suspected that this is related to the 
little-understood Pr ordering transition at 17K, but the connection is not clear. The 
parent Mott insulator PrB^CusOe has a similar phase transition at 11K (2D]- 

2. Pr systems 

In most compounds, Pr is trivalent. The existence of the dioxide PrC>2 shows that it can 
also be tetravalent though this is rare. The question of the valence state becomes 
crucial whenever the nature, or the driving force, of Pr ordering is in doubt. For Pr- 
filled skutterudites: the formation of a heavy Fermi sea, as seen in PrFe4Pi2 and in 
PrOs4Sbi2, is most likely to involve the hybridization of Pr /-states with conduction 
band states, and consequently deviation from strictly integral valence. 

These basic doubts notwithstanding, we propose to study the interplay of dipolar 
and quadrupolar interactions in localized /-electron models, and apply the results to 
PrFe4Pi2 and PrB^CusCV We assume the integer valence state Pr 3+ , and a Hund's 
rule ground state with J = 4. We note that a similar approach was successful in the 
case of CeB 6 pQ, though the system is known to possess the features of a dense Kondo 
system. We envisage a similar study of Pr systems. 

Pr 3+ is a non-Kramers ion. It follows that a unique ground state can be reached 
without breaking time-reversal invariance, either by having a singlet Ti crystal field 
ground state to begin with, or via the ordering of quadrupolar moments. 

2.1. The Ti-T 4 scenario for PrFe^Pu 

PrFe4Pi2 undergoes a second order transition at T tr = 6.51K, which is manifested 
in a susceptibility spike, and a A-shaped anomaly of the specific heat 0. Though 
at first thought to be antiferromagnetic, the ordered phase is now known to be 
antiferroquadrupolar. Upon switching on a magnetic field, the transition temperature 
decreases. The character of the transition switches to first order at the tricritical point 
Ttn ~ 5K, H tT[ ps 2T, and the ordered phase is completely suppressed at ~ 4T [2]. These 
parameters depend on the direction of the applied field. However, here we consider only 
the case of a field applied along the (100) direction. 
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The question whether the /-states are itinerant or localized can be posed again. 
PrFe4Pi2 is one of the few Pr-based heavy fermion systems [2|. Though having an 
even number of /-electrons per site makes the distinction between large and small 
Fermi surfaces difficult, we take the view that the large Fermi sea of heavy /-electrons 
competes with a state in which the conduction electrons build a small Fermi sea, and 
multipolar inter-shell interactions are switched on E3- This is supported by 

inelastic neutron scattering finding that the crystal field excitations become sharply 
defined when the temperature is lowered below T tr . It follows that we may use a localized 
/-electron model to describe the ordered phase and the vicinity of the phase boundary. 
We use the crystal field states of the octahedral group Oh- We note that the icosahedron 
of the twelve P atoms surrounding a Pr site gives a novel tetrahedral component to the 
crystal field potential pH], but we assume that using cubic symmetry labelling is an 
acceptable approximation. 

Since the local order parameter is one of the quadrupolar moments, it may look 
evident that the crystal field ground state should carry a (spontaneous) quadrupolar 
moment. However, an analysis of the high-T behavior of the magnetization curves shows 
that the level scheme (ground state)-(first excited state) may be either of the following 
three: Fi-Fi, Fx— Fs, or r 3 -r 4 [2j. Only the last, with the T 3 ground state, carries an 
unfrozen quadrupolar moment which can be ordered by switching on (arbitrarily small) 
intersite interactions. It was also shown that this choice is consistent with a symmetry 
analysis of the structural distortion accompanying the antiferro-quadrupolar ordering 
[HI] . This latter argument relies only on the assumption of the r 3 ground state, and 
does not refer to the nature of the first excited state. 

Here we show that the alternative Fi-Fi scheme is also capable to account for most 
of the observed static properties of PrFe 4 Pi 2 . We list the crystal field states 

T x = v /5/24(|4) + |-4)) + v ^7l2|0) 

r 4 ± = (l/4)(|3)±|-3) + v / 7[±|l) + |-l)]) 

r? = 072(1 - 4> - |4». (1) 

where we use the quadrupolar eigenstates Tf in the representation T 4 . Since the ground 
state carries no kind of moment, the ordered quadrupolar moment must be interaction- 
induced, implying that the interaction has to exceed a threshold value of the order of 
the ri-T 4 splitting. Taking it for granted that intersite interactions mix the T 4 states 
with Ti, the possible local order parameters are those appearing in the decomposition 

(r 4 © r 4 ) <g> (r 4 © r 4 ) = 2r x + r 3 + 3r 4 + r 5 . (2) 

The system is capable of either T 3 or T 5 type quadrupolar ordering, or any (dipolar or 
octupolar) kind of T 4 order. In principle, we have to allow for all the corresponding 
intersite interactions. However, we neglect the octupoles. We assume T 3 -type quadru- 
polar order and, recalling that the pair interaction has only tetragonal symmetry, and 
so 0\ and 0% are not equivalent, we arbitrarily keep only the 0\ = J% — Jy term. The 
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following mean-field-decoupled hamiltonian includes also a dipole-dipole coupling: 

W = 7~LcF + T^Zeeman + T^quad + ^dipole 




where we allowed for two-sublattice (A and B) order on the bcc lattice [z = 8). g — 4/5 
is the Lande factor, Q the quadrupolar, and I the dipolar coupling constant. We mention 
that a simplified interaction, with quadrupolar coupling only, was treated in a previous 
communication jTTj. However, the inclusion of 7^di P oie is important for getting a better 
fit to the experimental results, in particular the susceptibility. 




Figure 1. Left: The phase boundary in the H—T plane (H \\ x). Continuous line: 
second order, dashed line: first order transitions. Open circles: calculated, open 
rectangles: measured (after Middle and right: the inverse susceptibility l/x> 

and the susceptibility x, as a function of T. \ i s i n units of /j,B/Tesla/(Pr site). 

The diagonalization of (jSJ) in the basis (0) is straightforward. We show only the 
results. We have found that the (ferromagnetic) dipolar coupling I = 106mK, the 
antiferro-quadrupolar coupling Q = — 9.5mK, and the crystal field splitting A = 3K give 
good overall agreement with the observations (in case Q appears implausibly small, check 
that 0\ has big matrix elements). With these parameters, the only phase transition 
is the onset of antiferro-quadrupolar order at T tr = 6.5K. The transition is gradually 
suppressed by an external magnetic field (Fig. Q left). We locate the tricritical point 
at T tri pa 5K, H tvi pa 2T, in agreement with experiment. The same set of parameters 
gives also the observed Weiss temperature O = 3.6K from the intercept in the inverse 
susceptibility (Fig. [H middle). 

We note that our estimate of the crystal field splitting A pa 3K is substantially 
smaller than the one quoted in [2]; considering the scale of T tr , we have a quasi- 
quadruplet. However, this is not necessarily unreasonable: a calculation indicates a 
region where the gaps become very small (H. Harima, private communication). 

To conclude our review of the results from the r^-^ scheme!" , we show the 

| We did not make a comparable study of the full range of models, so we could not claim that the 
present model gives best agreement with experiments. Rather, we wished to demonstrate that the 
T1-F4 scheme is not obviously excluded, in spite of the absence of ground state quadrupolar moments. 
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Figure 2. Magnetization curves for the Ti-^ model, below (left) and above (middle) 
of the tricritical temperature, and in the disordered phase (right). Parameter values 
are the same as in Fig. ^ 



magnetization curves for three representative temperatures (Fig. In the region of 
first order transitions below the tricritical temperature, one finds a sharp metamagnetic 
transition which becomes smoother in the regime of continuous transitions. These results 
show a good overall resemblance to the measured curves 0. 



2.2. The role of quadrupolar coupling in PrBaaCu^Oe 

All the (RE)Ba2Cu307 systems, with the exception of RE=Pr, are Tc ~ 90K 
superconductors. Whenever the RE 3+ ion has a magnetic moment, they also undergo 
RE antiferromagnetic ordering at Tn = 1 — 2K. Analogous results hold for the parent 
Mott insulators (RE)Ba 2 Cu 3 6 , where most T N follow de Gennes scaling j2Zj, the peak 
belonging to RE=Gd with Tn ~ 2.2K |2Sj- This is well understood on the basis of the 
-^Gd-Gd = 156mK meaured by ESR [23]. We would expect that for RE=Pr, 

_ {(g-l?J(J+l)} Pr 



r N (Pr)| 



deGennes 



•T tr (Gd) » O.IK 



{(<?-i) 2 J(j + i)} Gd 

Instead, magnetic ordering and, at the same time, tetragonal-to-orthorombic distortion, 
is observed at T N « 17K for PrBa 2 Cu 3 07, and T N « 11K for PrBa 2 Cu 3 6 (201, exceeding 
the de Gennes estimate by two orders of magnitude. There must be a mechanism for 
producing substantially stronger exchange, or a different kind of ordering, or both. 

ESR on Gd:PrBa 2 Cu 3 6 gives the intersite exchange Ip r _Gd = — 140mK |2S|. De 
Gennes scaling would lead us to expect this sign, but also an absolute value of about a 
factor of 4 smaller. The greater spatial extent of Pr 4/ orbitals is likely to explain 
the enhancement of /; correspondingly, we would expect that /p r _p r is even more 
enhanced, but would probably still fall short of accounting for the observed Tn. Here 
we suggest that the presence of quadrupolar interactions may explain several features 
of the magnetic behavior of PrBa 2 Cu 3 06, and also an additional enhancement of Tn. 

All previous works agree that in the level scheme of Pr 3+ ions, a low-lying quasi- 
triplet consisting of the tetragonal doublet 



a| ±3) -p\ =f 1) 



(4) 
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and the singlet 



(|2> - I - 2)) 



(5) 



is well separated from the remaining six J = 4 states, and therefore suffices for 
modelling (as far as the Pr sites are concerned) all low-T phenomena ED] - With 
a = J7/8 ~ 0.9354, (J3J) and (J3J) would constitute the cubic T 5 triplet, but our fits 
yield a ~ 0.943, indicating a slight admixture from the doublet derived from T 4 . Let 
us observe that the doublet carries both J z dipole, and 0\ (or alternatively O xy ) type 
quadrupolar moment. If the doublet is the crystal field ground state, the degeneracy 
can be resolved either by magnetic, or by quadrupolar, ordering. 

Using (jlJ)-(jSJ), and following the procedure of j2l|, we fitted the high-T 
magnetization curves, and found that the dipole-dipole coupling has a strong planar 
anisotropy, and there is also a substantial ferroquadrupolar coupling of the 0\ moments 
[2E]. The latter finding may lead us to ask whether the observed transition is perhaps 
purely of quadrupolar nature; however, neutron scattering shows that T N is indeed a Neel 
temperature, with the T < Tn ordered moments strongly tilted out of the tetragonal 
c-direction |25| . Within the crystal field model, it is rather mysterious why the system 
does not take advantage of the permanent c-axis moments of the doublet (jlj), and chooses 
instead afe-plane moments which have to be interaction-induced. 

The mean field hamiltonian acting on the quasi-triplet is similar to © 



where we allow for two-sublattice antiferromagnetism, but assume ferroquadrupolar 
coupling. In view of the bilayer structure of PrE^CusOg, we take z = 4. From our 
susceptibility fits, we estimate A w 20K. 

We note that a closely related three-state problem was considered by Lfbero and Cox 
|2T] . They assume hexagonal symmetry, so their choice of the singlet and the doublet 
is quite different from ours, but nevertheless, when worked out fully, the results of their 
mean field theory should basically correspond to what we find. However, we calculate 
different quantities (such as the susceptibilities), and emphasize different aspects. 

Fig. El (left) shows a representative phase diagram. At Q = 0, the system gains 
energy from antiferromagnetic ordering, with the moments lying in the ±x-direction. 
However, it is easy to convince ourselves that at the same time, the system develops 
quadrupolar order: J x connects T° to Yf — Y^ (which is one of the quadrupolar 
eigenstates), but not to Y^ + Y^ (which would be the other). Thus, for zero or weak Q, 
the primary order parameter is (J x ), while is a secondary order parameter (Fig. El 
middle). At sufficiently strong Q, both interactions are important, and the two orders 
appear at a first order transition (Fig. El right). The regime of first-order transitions is 
bounded by a lower, and an upper, tricritical point (big black dots in Fig. El left). Above 
a A-dependent threshold value of Q/I, pure quadrupolar order sets in at a higher critical 
temperature. From this point on, the quadrupolar splitting is added to the crystal field 



"H — H-CF + ^Zeeman + T^quad + ^dipole 






Figure 3. Left: Transition temperatures as a function of the quadrupolar coupling 
Q[K], for / = 0.2K and A = 2K. Circles: calculated points; the dashed (first order) 
and continuous (second order) lines are drawn for convenience. Filled circles: tricritical 
points. Middle and right: magnetic (m) and quadrupolar (q) order parameters for 
Q = and Q = 0.05, resp. 



splitting, eventually suppressing magnetism. Note, however, that for a wide range of 
Q/I, magnetism is assisted by quadrupolar interactions, in the sense that the magnetic 
transition temperature increases with Q/I. This holds also for part of the regime where 
the two transitions are distinct. 




Figure 4. The T-dependence of linear (left) and nonlinear susceptibility (right) for 
purely ferroquadrupolar ordering. X3 diverges at the onset of quadrupolar order. 

A remarkable aspect of Pr ordering is its sensitivity to magnetic field, particularly 
for H _L (001) [221I2S1- Quadrupolar interactions are known to have a strong signature 
in non-linear magnetic response [21]. We show in Fig. 01 the results for I = 0, Q = 0.2, 
A = 2, for the field in the ab plane. The upward curvature of \, an d the divergence 
of X3 at the transition are reminiscent of the observed behavior. Fig. |S] shows the In- 
dependence of the magnetization for several fields H || (100). It is apparent that the 
sharp H = anomaly is quickly smeared out in higher fields, similar to the observed 
behavior of Prl^CusCV 

To conclude, we have examined models where the local Hilbert space is spanned 
by the states from the lowest two crystal field levels. The near-degeneracy of two 
levels has interesting consequences when the splitting is comparable to both intersite 
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Figure 5. Magnetization (/xe/Pr) versus temperature for the H | (1,0,0) magnetic 
fields 0.01T {left), 0.1T (center) and 0.4T (right). 

intractions, and laboratory magnetic fields. We found that the induced quadrupolar 
moment scenario gives a good understanding of the static properties of PrFe 4 Pi 2 , and 
that the inclusion of quadrupolar interaction is helpful in understanding the non-linear 
magnetic behavior of PrBa 2 Cu 3 6 . 
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